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Abstract 

The two-dimensional Ashkin- Teller model provides the simplest example of a statistical sys- 
tem exhibiting a line of critical points along which the critical exponents vary continously. 
The scaling limit of both the paramagnetic and ferromagnetic phases separated by the critical 
line are described by the sine-Gordon quantum field theory in a given range of its dimen- 
sionless coupling. After computing the relevant matrix elements of the order and disorder 
operators in this integrable field theory, we determine the universal amplitude ratios along 
the critical line within the two-particle approximation in the form factor approach. 



1 Introduction 



The quantitative description of the universality classes of critical behaviour is the main goal of 
quantum field theory when applied to statistical mechanics. For a statistical system possessing 
an isolated point of continous phase transition the canonical characterisation of the scaling 
behaviour of a thermodynamic quantity, say the susceptibility x, takes the form 



where t ~ J c — J measures the distance from the critical temperature 1/J C . Direct measures or 
numerical simulations of the system at different temperatures will then allow the determination 
of the critical exponent 7 and of the critical amplitudes r±. Contrary to the exponent, the 
amplitudes are not universal, but their ratio T + /r_ is pp. 

In principle, the universal quantities can be computed from the quantum field theory encod- 
ing the fundamental symmetries of the system. The critical exponents are yielded by the massless 
(conformal) field theory describing the critical point; for the amplitude ratios one needs instead 
working with the massive (the mass being an increasing function of \t\) field theory accounting 
for the deviations from criticality. 

Consider now the case in which the physical system exhibits a manifold of second order 
phase transition points on which at least some of the critical exponents vary continously. To be 
specific, we will refer to the simplest case in which the manifold is a line. Few remarks are in 
order about the comparison of measurements with theoretical preditions in this situation. The 
off-critical Hamiltonian of the system contains now two parameters J' and J" such that the 
critical line corresponds to a curve in the J'-J" plane. The field theory describing the scaling 
region close to the critical line will depend on a coupling (3 (which has nothing to do with inverse 
temperature) parameterising a line of fixed points of the renormalisation group, as well as on 
a mass scale measuring the distance from criticality. The correspondence between the values 
of (3 and the points of the critical line in the J'-J" plane is non-universal, i.e. it depends on 
the microscopic details of the system like the lattice structure. For each system in the given 
universality class, however, this correspondence can be determined comparing measurements and 
field theoretical predictions for at least one of the universal quantities (say a critical exponent) 
which vary continously along the line. Once this has been done, field theory yields predictions 
for all other universal quantities in much the same way as in the case of isolated critical points. 

For a given 0, the field theory describes a renormalisation group trajectory flowing away 
from the point that has been selected on the critical line. Calling t the coordinate along this 
trajectory, equations like (jl.lj) now define /3-dependent critical exponents and amplitudes. The 
trouble is that, in the generic case, it is not possible to locate the image of the given trajectory 
on the J'-J" plane. This means that the path along which the limit has to be taken in 
in not known and this kind of definitions are practically useless for measuring exponents and 
amplitudes in presence of continously varying critical behaviour. 

While the critical exponents can be measured in other ways (e.g. finite size scaling at 
criticality), the amplitudes appear essentially out of reach in the present case. For the purpose 
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of comparison with the field theory predictions, however, we are not interested in the amplitudes 
themselves, but rather in their universal combinations. If a duality transformation relating 
points on opposite sides of the critical line on the J'— J" plane is available, an amplitude ratio 
like r + /r_ can in principle be measured without actually measuring the amplitudes. Its value 
at a given point Pq along the critical line can be determined taking the ratio of the susceptibility 
measured at two dual points close to Pq, the theoretical error going to zero with the distance 
from Pq on the J'-J" plane. 

This paper deals with the field theoretical determination of universal ratios for the simplest 
class of continously varying critical behaviour. Critical phenomena in two dimensions are char- 
acterised by a number c, called central charge 2 , which increases with the number of degrees 
of freedom of the system. Since for c < 1 critical exponents are only allowed to take discrete 
values 0, the first possibility for continously varying exponents opens at c = 1. This is the 
central charge of a free massless boson (Gaussian model), which indeed possesses a continous 
one-parameter family of scaling operators. Different statistical models exhibiting a line of crit- 
ical points with continously varying exponents renormalise onto the Gaussian model at large 
distances. These include the Ashkin- Teller model [3] and the 8- vertex model which are 
related by a duality transformation j^j. Their precise relation at criticality with the Gaussian 
model has been determined in the past exploiting the exact Baxter solution of the critical 
8-vertex model. 

The Ashkin- Teller model is not solved on the lattice away from criticality 1 and the universal 
ratios cannot be computed exactly apart for the special case in which the model reduces to 
two decoupled Ising models. On the field theoretical side, while the relation of the scaling 
limit with sine-Gordon type deformations of the Gaussian model has been clear for longtime 
|S], the quantitative study has been prevented by the need of non-perturbative methods. Here 
we exploit the integrability of the sine-Gordon model to compute the universal ratios along the 
Ashkin- Teller critical line in the two-particle approximation within the form factor approach. 
The form factors of the sine-Gordon model have been and continue to be the subject of intensive 
study (see [§1 1101 ITT1 112j among other references). The description of the Ashkin- Teller model, 
however, relies essentially on the control of the order and disorder operators a and /i which are 
not among those considered in these works. We compute all the one- and two-particle matrix 
elements needed for our purposes within the framework based on the properties of mutual 
locality between particles and operators. Duality is used to describe through the sine-Gordon 
field theory both the paramagnetic and ferromagnetic phases on the two sides of the critical line. 

The paper is organised as follows. We review the phase diagram of the Ashkin- Teller model 
in section 2 and its field theoretical description in section 3. Section 4 deals with the scattering 
theory for the scaling limit around the critical line while section 5 is devoted to form factors. 
Correlation functions are discussed in sections 6 and universal ratios are computed in section 7 
before few concluding remarks. 

1 More precisely, away from its self-dual line. 



2 



2 The isotropic Ashkin- Teller model on the square lattice 



The Ashkin- Teller model (3J I13| describes two Ising models coupled by a four-spin interaction. 
The case with the two Ising models having the same temperature 1/J is called 'isotropic' and 
corresponds to the Hamiltonian 

Hat = ~ J2{J[ai(x)ai(y) + a 2 {x)a 2 {y)) + Ji(Ji(x)ai(y)a 2 (x)a 2 (y)} , (2.1) 

where &i{x) and cr 2 (x) are the two Ising spins at site x and the sum runs over nearest-neighbour 
pairs {xy). Each of the transformations o\ — > — a\, a 2 — > — a 2 and o\ <-» a 2 leaves the Hamilto- 
nian invariant. 

We are interested in the square lattice model for which the transformation J — > — J amounts 
to reversing the spins o\ and <r 2 on one sublattice. In this case the phase diagram is symmetric 
under reflection about the J4 axis, under which ferromagnetic ordering in o\ and a 2 becomes 
antiferromagnetic, and vice versa. With this remark in mind we will only refer to the case J > 
in the following. 

Obviously, the model possesses a critical point in the Ising universality class at J = J* = 
^ln(l + \/2) along the decoupling line J4 = (the point marked D in Fig. 1). The model 
becomes invariant under permutations of the four states (a±, a 2 ) along the line J = J4 which is 
then the 4-state Potts subspace. Hence, a Potts critical point is located at J = J4 = \ In 3 (point 
P in Fig. 1). The antiferromagnetic 4-state Potts model on the square lattice is non-critical |14j . 
so that no other critical point for J > is implied by the analysis of the Potts subspace. 

When J = 0, the Hamiltonian (|2.1j) describes an Ising model in the variable <Ji<j 2 . Thus the 
points If and Iaf in Fig. 1 located at J4 = ±J* along this line are a ferromagnetic and an 
antiferromagnetic Ising critical point, respectively. 

As J4 — > +00 the energies a\{x)a\{y) and cr 2 (x)a 2 (y) take the same value, their product 
being forced to 1. Thus in this limit the Hamiltonian (|2.1j) reduces to a single Ising model with 
coupling 2 J. A ferromagnetic Ising critical point is then located at (J, J4) = ( J*/2, +00). 

The model admits a duality transformation which maps the Boltzmann weights 

e 2J+J 4 +J 
g— J4+J0 
e -2J+J 4 +J 



(Jo is an arbitrary constant) into the new 

tbi = 
uj 2 = 

U) 3 = 



UJl = 
U>2 = 
UJ 3 = 



ones 
1 

-(uJi + 2uj 2 + uj 3 ) 

-(ui-u 3 ) 

1 . . 
-(uji - 2uj 2 +uj 3 ) . 
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Figure 1: Phase diagram of the isotropic Ashkin-Teller model on the square lattice. The self- 
dual curve is divided into a critical line with continously varying exponents (continous line) 
and a non-critical part (dotted line). The dashed curves are critical lines with Ising critical 
exponents. The dash-dotted line is the 4-state Potts model subspace. Four different phases are 
labelled by the roman numerals. 
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For J > the weights tOj are positive (and the corresponding couplings J, J4 and Jo real) 
provided 

cosh 2 J > e" 2j4 . (2.2) 

Within the region selected by this condition duality relates points located on opposite sides of 
the self-dual line determined by the equation 

sinh2J = e - 2j4 (2.3) 

(the line A-D-P-B in Fig. 1). 

A duality transformation on the a\ spins only maps the square lattice Ashkin-Teller model 
onto a staggered 8- vertex model [S]. In the isotropic case (|2.1|) the staggering disappears along 
the self-dual line (|2.3|) . Since the unstaggered 8- vertex model is exactly solvable, Baxter [3] 
was able to show that the self-dual line is critical for J4 < ^ln3 (curve A-D-P in Fig. 1) and 
non-critical for J4 > |ln3 (curve P-B in Fig. 1). The critical exponents vary continously along 
the critical portion of the line. We mention that the value J4 = J*/2 selects on the critical line 
the Fateev-Zamolodchikov Z4-parafermionic critical point |15j . 

No exact lattice results are avalaible for the model away from the self-dual line. The complete 
phase diagram, however, has been obtained through a variety of approximate methods (see |16j). 
The critical line with continously varying critical exponents bifurcates at the Potts point P into 
two critical lines, dual of each other, ending into the previously mentioned Ising critical points 
(J, J4) = ( J*/2, +00) and Ip. Another critical line originates from the antiferromagnetic Ising 
critical point at Iaf an d points towards J4 — * —00. The exact location of these three critical 
lines (dashed in Fig. 1) is unknown. The critical exponents are expected to stay within the Ising 
universality class along them. 

The four critical lines are the boundaries of four different regions in the phase diagram of 
Fig. 1. In I the system is ferromagnetically ordered with (01}, (02) and (01 02) an different from 
zero; phase II is the disordered one in which all these order parameters vanish; phase III exhibits 
partial ferromagnetic ordering, with (a\) = (02) = but (01 02) 7^ 0; phase IV is similar to III 
but with 0102 ordered antiferromagnetically. 

3 Scaling limit and field theory 

From the field theoretical point of view the critical line with continously varying exponents of 
the Ashkin-Teller model must correspond to a line of fixed points of the renormalisation group 
generated by a marginal operator. The latter can be easily identified looking at the point D 
on the line where the model reduces to two non-interacting critical Ising models. The four-spin 
term in (|2.1jl gives in the continuum limit the product Ei£2 of the energy densities of the two 
Ising models. Since the scaling dimension of the energy operator in the Ising model is 1, £162 
is indeed marginal. The existence of a line of fixed points implies that this operator remains 
strictly marginal even away from the decoupling point and is responsible for the deviation of the 
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critical exponents from the Ising values. Accordingly, we write the action for the scaling limit 
around the line of fixed points as 

•Ashling = -A® + A® + T J d 2 X (e X (x) + E 2 (x) ) + p Jd 2 X£ 1 E 2 (x), (3.1) 

where A®, j = 1,2 denotes the fixed point Hamiltonian of the j-th Ising model. The couplings r 
and p are non- universal functions of the lattice couplings J and J4. They determine the distance 
from criticality and the coordinate along the fixed line, respectively. 

The Hamiltonian (j3.1j) implies that the central charge along the line of fixed points is twice 
that of the Ising model, namely c = 1. This agrees with the fact that the 8- vertex model along 
its critical line reduces to the 6- vertex model, and the latter renormalises al large distances on a 
free massless boson, namely the Gaussian model with central charge equal to 1. It is well known 
|171 118j that a Gaussian fixed point in two dimensions also admits a description in terms of a 
Dirac fermion ip = ipi + iip2- The neutral components ipj are the Majorana fermions associated 
to the j-th. Ising model. The energy operator Ej is bilinear in the fermion (sj = "tpfipj , with 
tjjj denoting the two components of the spinor), while the spin operator crj and the disorder 
operator p,j are non-local with respect to it. 

The study of the Ashkin- Teller critical line in terms of the Gaussian model was performed 
in [7]. At the Gaussian fixed point, the boson field can be decomposed into its holomorphic 
and antiholomorphic parts as <f(x) = <fi(z) + 4>(z), where we introduced the complex coordinates 
z = x\ + 1x2 and z = x\ — 1x2- The scaling operators of the theory are the vertex operators 

V PtP {x) = jiP^+vim i (3.2) 

with conformal dimensions (A, A) = (p 2 / '8tt , p 2 / '8tt) , scaling dimension X = A + A and spin 
s = A — A. They satisfy the gaussian operator product expansion 

( x )^P2,p2 (0) = z PlP2 ^ zPiPz/ Vp 1 +p 2j p 1 +p 2 (0) + . . . . (3-3) 

We see from this relation that taking V Pli p 1 (x) around V P2t p 2 (0) by sending z — > ze 2t7T and 
z — ► ze~ 2l7T produces a phase factor e 2i7r71 ' 2 , where 

7i,2 = ^- ~ P1V2) (3.4) 

is called index of mutual locality. If 7^2 is an integer the correlator {Vp lt p 1 (x)V P2t p 2 (0)) is single 
valued and the two operators are said to be mutually local. Since 71 1 = 2s, the operators which 
are local with respect to themselves (the only ones we are interested in here) must have integer 
or half integer spin. 

Without loss of generality, we call cos j3tp the (most relevant component of the) energy 
operator £ = e\ + £2 which drives the system away from criticality. Here j3 is the parameter 
which accounts for the continously varying exponents and then parameterises the critical line; 
(3 2 is equal to 4-7T at the decoupling point where the scaling dimension of E must be equal to 1. 
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All the operators of interest for the description of the lattice model must be local with respect 
to the energy operator. This locality requirement selects the operators V PtP with p — p = Airm/ (3, 
m integer, namely 

V P)P {x) = e^W (3.5) 

and 



U n , m (x) = Vn0, a- &_&Jx) = j^^+T™^} , n = 2 S = 0, ±1, . . . , m = ±1, . . . 

2m (3 ' 2m (3 

(3.6) 

Here we introduced the 'dual' boson field <p which is 4>(z) — (j>(z) at criticality and satisfy the 
relation 

■ 9<f d<p 

= £ ab^— ■ ( 3 - 7 ) 

ox a dx b 

The operators e iptp and ?7o,m are scalars (s = 0) and have scaling dimensions X p = p 2 /Air 
and Xq )JTI = 7rm 2 /P 2 , respectively. The lowest operators with |s| = 1/2, i.e. £/±i,i and U±\-\ 
with conformal dimensions A„ jm given by 

A ±u = A ±w = -(£- ±2 + _j (3 . 8) 
A ±M = A ±1 ,_ 1 = ig T 2 + ^), (3.9) 

form the Dirac spinors 

\<l>- J \<K+ Hq ) \U-i,i) 

\ </>'-) \ *: - « I v / 

The Ashkin-Teller operators and their bosonic form [7j are listed in the first two columns 
of Table 1; the scaling dimensions are given in the third column. The rest of the table specifies 
the behaviour of the operators under the following tranformations: 

Spin reversal . We refer to the invariance of the model under the reversal of all o\ or all 02 
spins as Z2 X Z2 symmetry. The transformation 

J<P - «-j<P 

of the bosonic fields is found to correspond to the reversal of the u\ spins. The transformation 

tp^-tp 
<p^-<p 

corresponds instead to the reversal of the 02 spins. 
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Bosonic form 
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Table 1: Ashkin- Teller operators with their bosonic counterparts, scaling dimensions and sym- 
metry properties. When a symmetry transformation sends an operator into ± itself, only ± is 
indicated. 



Exchange . The exchange of the two Ising copies (1 <-> 2) is implemented in the bosonic 
language by the transformation 

2"7T _ JT 27T . 

Semi-duality . We call semi-duality the transformation Dj which interchanges Oj and jij. We 
saw in the previous section that a transformation of this kind relates the Ashkin-Teller model 
to the 8-vertex model. In the bosonic language D\ corresponds to the exchange 

(3-12) 

and L>2 to 

2^-2"^^ (3-13) 

Clearly Dj = 1. Notice that D\ is not quite the duality transformation of the first Ising model. 
In particular, it changes the sign of £2 rather than £\. A similar observation applies to D2. 
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Duality . This is the full duality transformation of the Ashkin-Teller model and is obtained 
composing D\ and D<i- Due to sign factors there are two possibilities, D + = D1D2 and D- = 
D2D1. D± correspond to the bosonic transformations 



(3(f —>■ (3if ± 7T 

2-7T 2tt it 



Notice that D + D_ = D_D + = 1, 1. 



3.1 The self-dual line 



The field theory describing the self-dual line (|2.3|) at large distances has to be invariant under 
all the above transformations. This requirement selects the action 



A sd = J d 2 x ^ (d a ip) 2 - {dn cos 2n/3<p + g n cos 8n^<£ 



(3.14) 



where the g n 's and g n 's are, like (3, non-universal functions of J4. This theory is critical (massless) 
as long as none of the operators perturbing the Gaussian term becomes relevant. Since cos 2/3<p 
becomes marginal at (3 2 = 2tt and cos87rc/3//3 at f3 2 = 8tt, the critical range is 2tt < j3 2 < 8tt. 

The 4-state Potts critical point (P in Fig. 1) is the right end point of the Ashkin-Teller 
critical line. At this point the Ising variables a%, 02 and <J\o% play a completely symmetric role 
and must have the same scaling dimension 1/8. Hence, it follows from Table 1 that the Potts 
point corresponds to (3 2 = 2ir. 

The relation between (3 and J4 along the critical line for the case of the square lattice model 
with nearest neighbour interactions can be obtained comparing the energy scaling dimension 
predicted by the Gaussian theory with that coming from the lattice solution. It reads [7j 

4vr i 2 . ( tanh2J 4 \ 

■ arcsm +„„i.ot T > ( 3 - 15 ) 



(3 2 ir Vtanh2J 4 -l 

from which we see that the limit J4 — > —00 corresponds to (3 2 = 6tt. Thus the Ashkin-Teller 
critical line with continously varying exponents spans only a portion of the critical region of the 
theory (|3.14j) . namely the range 

2vr < (3 2 < 6vr . (3.16) 

For (3 2 < 2tt at least one of the operators cos2n(3ip is relevant. The theory (J3.14|) is massive 
and corresponds to the non-critical part of the self-dual line. The relation (|3.15|) does not hold 
in this region. 



3.2 Breaking duality 

The field theory describing the model on the two sides of the self-dual line is obtained adding 
to the self-dual action (|3.14|) the operators cos(2n — l)f3ip, which preserve the spin reversal and 
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exchange symmetries but are odd under duality 

Aat = A sd -^Tn d 2 x cos(2n - l)/3tp , (3.17) 
n •* 

where the r n 's are functions of the lattice couplings J and J4. 

To describe the scaling regions around the critical line with continously varying exponents 
we keep the only operator which is relevant in the range (|3.16|) . Thus we are left with the 
euclidean sine-Gordon action (we set t\ = t) 

A SG = J d 2 x Q {d a v? - t cos (3^j (3.18) 

which is the bosonic version of (|3.1|) . 

For r > this action describes the scaling region of the paramagnetic phase II in which the 
vacuum |0) located at (p = is invariant under spin reversal and exchange symmetry and we 
have 

(OltTi 1 0) = <0|<7 2 |0> = 
<0|o-i<t 2 |0) = (0|sin^|0> = . 

For r < the two vacua |0±) located at 0(p = ±tt are the image of the high-temperature 
vacuum |0) through the duality transformations D±. The vacua |0+) and |0_) are interchanged 
by the spin reversal and exchange transformations, so that the internal symmetries of the model 
are spontaneously broken (ferromagnetic phase I). We have 

(0+M0+) = (0+M0+) = ±M a 

<0_|(7i|0_) = — <0_|cr a |0_> = ±M a 
(0 ± |<7ia 2 |0±) = (0±| sin ^|0±) = ±M V 

where M a and M-p are positive. The last equation follows from the relations 

(0+|e iQ ^|0+) = (0_|e~^|0_) 
(0 ± |e iQV |0 ± ) =e ±2i7ra//3 (0±|e-^|0±) . 

3.3 Bifurcation at the Potts critical point 

For P 2 < 2tt more operators become relevant in the action (|3.17|) . Neglecting all irrelevant terms 
in the range 8tt/9 < (3 2 < 2tt leads to the double sine-Gordon action (g± = g) [5| 

Adsg = Jd 2 x Q (d aV ) 2 - t cos 0<p - g cos 2/3^J . (3.19) 

This quantum field theory has been analysed in Ref. ^Hj. The mechanism through which it 
accounts for the bifurcation of the critical line at the Potts critical point is easily understood 
already at the classical level. To see this let us fix a value of j5 in the considered range and 
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Figure 2: Schematic phase diagram of the double sine-Gordon quantum field theory (|3.19j) for 
fixed (3 2 < 2tt. Two massless trajectories starting from the Gaussian fixed point at the origin flow 
towards infrared Ising fixed points and divide the plane into three different phases. The points 
at ±r c on these trajectories correspond to dual points along the Ising critical lines bifurcating 
from point P in Fig. 1. 



treat for a moment both r and g as free parameters. For g = and fixed r > the vacuum 
of the theory (|3.19j) (i.e. the minimum of the potential invariant under the symmetries of the 
paramagnetic phase) is located at ip = 0. It stays there as we decrease g down to a critical value 
g c , classically equal to — r/4, where the quadratic term of the potential vanishes. For g < g c 
we have a maximum at ip = and two new minima located symmetrically with respect to it. 
Thus an Ising phase transition occurred at g c . Similar considerations can be repeated for r < 
(starting with one of the vacua at ip = ±7r for g = 0) and the resulting picture can be confirmed 
at the quantum level |19j . In summary, the lower r-g half plane is divided into three regions 
(Fig. 2) by two massless trajectories (corresponding classically to g = ±t/4 but whose precise 
location is unknown in the quantum theory) along which a flow from the Gaussian fixed point 
at the origin to infrared Ising fixed points takes place. Region II is a paramagnetic phase where 
the vacuum is located at ip = 0; region I is a ferromagnetic phase dual to II where ((3(p) equals 
7r or — 7r; in region III {(3ip) interpolates smoothly from to ±7r taking the value ±7r/2 at r = 0. 
In this latter phase we have (cxi<72) 7^ 0; it is less clear to us how to argue that (<7i) = (02) = 0. 

In order to make contact with the phase diagram of Fig. 1 we need to recall that g is not a 
free parameter at fixed j3. Indeed, both and g are determined by the value that J4 takes along 
the self-dual line. It is the value g{(3 2 ) that determines the distance r c of the Ising transition 
points from the given point along the self-dual line (see Fig. 2). The phase diagram requires 
that g{0 2 ) vanishes at (3 2 = 2ir and then decreases with (3 2 . The result in the coupling space of 
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Figure 3: Like in Fig. 2 but with varying /? 2 . The oriented lines indicate the flow towards larger 
distances along the critical surfaces. The thick lines are determined by the intersection of these 
critical surfaces with the image of the J-J4 plane of Fig. 1. 

the theory (j3,19|) is shown in Fig. 3. For 7r/2 < j3 2 < 2ir the non-critical part of the self-dual 
line is described by the action (|3.19|) with r = 0. As the Potts critical point is approached from 
below along the line {fi 2 — ► 2ir~ , g — > 0~) the operator cos2/3(/? becomes marginally relevant 
implying an exponential rather than power law divergence of the correlation length. The same 
massive field theory describes the q — > 4 + limit of the (/-state Potts model at T = T c [2"T)] . 

4 Scattering theory 

According to the discussion of the previous section the scaling limit of the Ashkin-Teller model 
around the critical line with continously varying exponents is described by the sine-Gordon 
theory (|3.18|) . This quantum field theory is integrable and the associated elastic and factorised 
scattering matrix is exactly known ([2^ an d references therein). 

The elementary excitations are the soliton A + and antisoliton A_ which interpolate between 
adjacent vacua of the periodic pontential. While being topologic excitations of the bosonic action 
(j3.18|) they correspond to the fermions if) and ip* of the equivalent fermionic model (the massive 
Thirring model) |17| ll8j. Actually, the integer m in (|3.6|) measures precisely the topologic charge 
and all the operators with m = 1 (—1) create a soliton (antisoliton) when acting on the vacuum 
of the theory. Writing A± = {A\ ± iA?)l^f2, the operators ipi and criAt^ +1 ( mo d 2) are both 
suitable interpolating operators for the neutral component A{. To be definite, we will refer to 
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Figure 4: Simple pole diagram associated to Eq. (|4.7|) . 



the latter choice in the following (see Table 2). The scattering of the particles A\ and A 2 in the 
integrable theory is completely determined by the relations 2 



A i (9 1 )A j (9 2 )= S i-A 



02)M6 2 )A k ( 



fcj=l,2 



with the non-zero scattering amplitudes given by 21 



where 





= Sj±(9) = 


s(9) + s+(e) 
2 






S+(9) - S(9) 
2 






S(9) + S-(0) 
2 






S(0)-S-(0) 
2 



S+(9) 
S(P) = 



S(9) 



sinh ^(6> + m) 
sinh — in) 
cosh -^(9 + Z7r) 
'cosh^^-^) 5 ^ 



cxp 



x sinhf^cosh| 



0x 



sm ■ 



7T 



vr/? 2 



8vr-/3 2 



(4.1) 



(4.2) 



2 The on-shell energy and momentum of a particle of species a are parameterised as (p°,p) 
(m a cosh 9, m a sinh 8) . 
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For (3 2 < Air (attractive regime) the above amplitudes possess simple poles in the physical 
strip Im6 6 (0,7r) corresponding to bound states (breathers) B n . The poles are located at 
6 = i(ir — n£) (n even for S + (9) and odd for S-(6)) and determine the masses of the breathers 



m„ = 2M sin ■ 



1< n < 



(4.3) 



where M is the mass of the particles Aj and [x] denotes the integer part of x. The lightest 
breather B\ is the particle interpolated by the boson ip and then is odd under all the internal 
symmetries of the model. The breather B n can be seen as a bound state of n breathers B\, so 
that its symmetry properties are those summarised in Table 2. 



Particle 


Creating 
operator 


Z 2 x Z 2 


1 2 


A 1 


fl/X2 


- x + 


A 2 


A 2 




+ x - 




B 2 k+1 


sin/3(/9 


— X — 




B2k 


cos (3ip 


+ X + 


+ 



Table 2: Particles, intepolating operators and their symmetries. 

The scattering of the breathers with the elementary excitations and with themselves is com- 
pletely diagonal (initial and final states are identical). The corresponding amplitudes are 

n5S 1,/2 *| +( „_ 2 ^W, noM 
Sj,8,(«)=tu=s(»)=<( (1-4) 



^)nfc 2) % {n _ 2j)f (0), neven 



where 



min (m,n)-i 

tanh±(fl + i7ra) 

*a(0) = 7 r : — r • (4-6) 

tann ^ (ft — 27ra) 

The simple pole at = i(7r + n£)/2 in the amplitude SA 3 B n (&) corresponds to the appearance 
as a bound state in the channel AjB n of the particle A, + i(mod2) ( n odd) or (n even). The 
simple pole at = i(n + m)£/2 (0 = i[7r — |n — m|£/2]) in the amplitude SB m B n (9) signals 
that B n+m (-Bi n _ m i) appears as a bound state in the B m B n channel. A crossed channel pole 
located at m — 9 is associated to each of these direct channel poles. A direct channel pole at 
6 = iu c ab corresponding to the particle c appearing as bound state in the ab channel allows the 
determination of the three-particle coupling T c ab through the relation (see Fig. 4) 

S ah (6 ~ < b ) ~ i^^- (4.7) 
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Figure 5: Double pole diagram associated to Eq. (|4.8j) . 
where Tf = C a C b C c T c ab , T c ab = C a C c Y a bc with C Aj = 1, C Bn = (-1)™. We also note the property 

r i n A = (-i) i+i rA n A 

The second order poles do not correspond to bound states and can be explained in terms 
of multiscattering processes of the type shown in Fig. 5 |22[ I23| . In the vicinity of such a pole 
located at 9 = iip the scattering amplitude can be written as 

S ab (0 * Up) * r ^ r ^ rg _ d ^f ce(i7?) . (4.8) 

5 Form factors 

The knowledge of the scattering amplitudes allows the determination of the form factors 

^,..a„(#i, ■ • • A) = <O|$(O)|ai(0i) . . . a n A)) . (5.1) 

Many results are known about the form factors of the sine-Gordon model 1101 1111 112j . Here 
we will recall how the one- and two-particle form factors can be computed within the general 
approach. These matrix elements for the order and disorder operators er, and fij, which are 
crucial for the description of the Ashkin-Teller model, were first considered in |24j . 

We use for the two-particle form factors of a scalar operator $(x) the notation 

F%(e 1 -e 3 ) = (oMo)\a(o 1 )b(9 2 )), (5.2) 

the dependence on the rapidity difference being a consequence of Lorentz invariance. The matrix 
elements satisfy the equations El IZEj 

Fl(9) = S c a i(8)Fg(-9) (5.3) 
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me + 2i7r) = h, a F^(-9) (5.4) 
Res e=i7r if 6 (0) = i5 ab (1 - (5.5) 
Res e=iul F*(0) = ir c ab F? (5.6) 

where the factor l$^ a takes into account the mutual locality between $ and the particle a (see 
Table 3). If <fi a is the operator which interpolates the particle a, then = e 2171 " 7 *-'^ . For the 
operators having a bosonic expression, 7$ i9 i a is given by (|3.4[) . For <7j- and /ij the results of 
Table 3 follow from the properties la h a k = ^j,^ k = 1 an d ^o-.,-,/^ = (— l)* 5 ^ (notice that the even 
breathers appear in the AjAj channel while the odd ones appear in the A\A<i channel). 





Bosonic 
form 












("If 








("If 


8 = ex + e 2 


cos flip 


+1 


+ 1 


C = e\ - e 2 


COS -jg-(^ 


+1 


+ 1 


£+ = -0x^2 + "02^1 


47T ~ 

sm 


+1 


+ 1 




sin/3</9 


+1 


+ 1 


V = <J\<J2 


sm 


-1 


+ 1 


V* = 11X112 


COS ij(p 


-1 


+ 1 




COS 


+1 


+ 1 


A*lO"2 


27T ~ 

sm ^93 


+1 


+ 1 



Table 3: Factors of mutual locality between particles and operators entering Eqs. Q5.4JI and (|5.5|) . 



The form factors are further constrained by the asymptotic bound 26 

lim F%(6) ~ e^l l 



(5.7) 



(5. 



where X$ is the operator scaling dimension. The non-zero two-particle matrix elements on the 
elementary excitations A± and A2 in the disordered phase are then uniquely determined to be 

m(ji) F (6) 



2£o;(mt) sinh ^{9 — m) 

F (e) 



-1 



cj(2ivr - 9) 



cosh I 



F5.,.(0) = i Cl -. , . 
^£a,W = c 2 (-1) j FoW 



27T 



cosh ■ 



-ci 



cosh || ( 



Z7T I 



F o (0) 



(5.9) 

(5.10) 

(5.11) 
(5.12) 

(5.13) 



16 



7T 



F {9) 



ITT 



£ cosh jg(0 — in) 
F (9) 



£ x ' sinh^(6»-i7r) ' 
Here c\ and C2 are normalisation constants, (fx) = = {^2), and the functions 

dx sinh (l ~ I) x sin 2 |f 



w(0) = exp { 2 



J 



x s inh ^ sinh 2x 



^b(#) = — isinh- exp ^ j : 



dx sinh § (l " I) sin 2 



(i7T — 6)x 
2^ 



x sinh f| cosh | sinh x 



satisfy the equations 



sinh£(6> + i7r) 

u(9) = u>(-6), oj(8 + 2iTr) = - . * , M 9-2i7r) 

Sinn 2({9 — iTT) 

F (9) = S(9)F (-9) , F (# + 2ivr) = F o (-0) . 
For large values of \9\ they behave as 



oj{9) ~ exp 
Fq(9) ~ exp 



7T 



+ 1 



For £ < 7r/2, the analytic continuation 

sinh[(0 - i(7r - 2£))/4] sinh[(fl + i(ir - 2£))/4] 



u;(2ivr - 0) 



cos 2 [(^-2£)/4] 



sinh f 1 - f ) » „ (2ivr - 0) 
^ — sin - 



x exp < 2 / — 

'0 x s inh # sinh 2x 



is convergent for real rapidity values. 

The breather-breather form factors can be written in the form 



2vr 



FIbJ9) = PIbJ9) cosh 



0^(i-'*, B „)(-i)*"' m F™ n Bm {8) 

F>B n B m {9) 



Here 



min (m,n)-i 

J'=l 



The function 



T a (9)=exp{2 



dt 



cosh (a - \) t f i7r _ e \ t 
-\ sm 2 - — - — - 



t cosh I sinh t 



2vr 
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solves the equations 

T a {9) = -t Q (6)T a (-0) , T a {9 + 2m) = T a {-9) (5.26) 
and behaves asymptotically as 

T a (0) ~ exp(|0|/2) , \9\^oo. (5.27) 

In particular 

T o (0) = Ti{p) = -isinh^ . (5.28) 

The denominator DB n B m {&) in (|5.23j) accounts for the pole structure of the form factors and 
can be written as (see |26| ) 

min (m,n)-i 

D Bm B n (e)=V im+n)4 _{e)V 1 _ lm _ nl j_(9) ]J >V-n| + 2,)xW (5-29) 

J=l 

with 

_ ,. x cos ttq. — cosh 6 , „ ,„ „„. 

P Q « = —5 2^ » « # 1 5.30 

2 cos 4p 

Pi(0) = l (5.31) 

W«(0) =P a (0)Pi_ Q (0) . (5.32) 
The last ingredient of (|5.23|) are the polynomials 

N ab 

^W = Ec a ?(cosh0) fc (5.33) 

k=0 

whose total degree is constrained by the the asymptotic bound (|5.8j) and whose coefficients are 
determined by the residue equations (|5.5j) and (|5.6j) together 3 with the relation |2H] 

associated to the double poles (|4.8|) in the scattering amplitudes. In particular, one finds 
iV| iBi = and N% 2B2 = 1 for <D = M , £ , 7>*, and JV^ ljBa = 1. 
Concerning the channel AjB n , the property 

= (-i) Bf W*) ( 5 - 35 ) 

has to be taken into account. We find that the minus sign appearing for odd breathers (and 
implying ^a^b^i"^ = — ^b^^a^ 2> ) * s nee ded if solutions compatible with the asymptotic 
bound ()5.8|) are to be found for the operators Oi and Cj/Lt i+1 (modi)- Then one can write 

/ a\S n ,oddSi,u B FTV, (9) 
<B,W = Fl i B„W (cosh |) '"^y. (5-36) 



J See also l)6.rifl for the energy operator £(x) 
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where 

n— 1 



2 r(n-2j) 2iv 
j=0 



n& 1)/2 ^ + M)i(«). «o 

D A] B n {e) = v l+ ^m < { (5.38) 

P|Wnt 2)/2 ^ +M) iW, neven 

and P% B n {0) are polynomials in cosh# to be fixed through the conditions on the poles. In 
particular we have 

F ^ + u mod 2 ){e) = _ F ^ {mod2){0) = { _ iy a r (1+ g /7r)/2 (g) ( ^ 39) 



dd 



J J n/2if)n/2-K/fl f, J 

with a, 6 and c determined by the equations 

Res e=i(7r+s)/2 F^ Bi (0) = ir^ B F- (5.41) 
Res, = , (7r/2+o F^ B2 (0) = iT^Fg (5.42) 
Res e=i7r/2 F^ B2 (0) = iT% B T B B \ B F% Bi {i{* - 30/2) . (5.43) 

The last equation is the specialisation of (|5.34j) to the double pole appearing at 9 = in/2 in 
the AjB2k scattering amplitudes and related to the diagram of Fig. 5 with a = Aj, b = -B 2 fc, 
d = e = B}~ and c = Ai, i = j (i ^ j) for k even (odd). It can be checked that the matrix 
elements ()5.40(l determined in this way satisfy the asymptotic factorisation condition [2S] 

tefW)-*^ < 5 - 44 > 

with k = 1. The condition 



lim 



fixes the relative normalisation between order and disorder operators. 

6 Correlation functions 

Within the form factor approach, correlation functions are obtained through the spectral sum 

<*i(x)* 2 (o)> = E E / fr-f^ 

<O|$i(O)M0i) . . . a n (9 n ))(a n (9 n ) . . . a 1 (0i)|$ 2 (O)|O)e- B ™l a: l , (6.1) 

where 

n 

E n = E m a fc cos h #A; (6.2) 
k=l 
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Figure 6: Different approximations of the central charge coming from Eqs. and (|6.1|l . The squares 
indicate the contribution of the states AjAj only, the diamonds the inclusion of the states B-i and B\B\, 
the stars the inclusion of the states B4, B1B3 and B^B-i- The exact result is c = 1. From right to left, the 
dashed vertical lines correspond to the first four thresholds where the breather B n (n = 1, . . . ,4) enters 
the spectrum of asymptotic particles. 



c 



0.995 - 
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Figure 7: A detail of Fig. 6. 
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denotes the total energy of the n-particle asymptotic state. This large distance expansion can 
produce good approximations of the integrated correlators even when only few lightest states 
are included in the sum. A quantitative illustration of the convergence pattern is obtained using 
the spectral sum to compute exactly known quantities like the central charge c and the scaling 
dimensions X$ through the sum rules [23 EE] 



J d 2 x\x\ 2 (Q{x)e(0)) c (6.3) 



3 

c = - — 

Air 

where (• • -) c denotes connected correlators and O(x) is the trace of the stress-energy tensor. The 
latter is proportional to the energy operator £(x) and its form factors are normalised through 
the condition 

F®{m) = 27rm 2 a . (6.5) 
Moreover, for a ^ b, P^{0) factorises a term j^Hl 

coshe+ m } + m2b . (6.6) 



The asymptotic condition [2"H] 



2m a m b 



lim F§ B (9) = (6.7) 

e^oo BnBn v ; (0) v ; 



can be used to determine the expectation value 



(9) = -irM 2 tan | , (6.8) 



a result which coincides with that known from the thermodynamic Bethe ansatz (see |29|). 

Figures 6 and 7 show the first few approximations provided by the insertion of a truncated 
spectral sum into the formula ()6.3|) for the central charge. For this quantity as for the scaling 
dimensions computed through l|fi.4j) . the states AjAj yield the exact result at the free fermion 
point f3 2 = Air, as well as the state B\Bi gives the exact result at the free boson point /3 2 = 0. 
Away from these free points the convergence of the series is extremely rapid due to the factor 
\x\ 2 in (|6,3|) which suppresses the contribution to the integral of the short distances, namely the 
region where the truncated spectral series fails to reproduce the exact correlator. Hence, in the 
repulsive region 4-7T < j3 2 < 8tt the AjAj contribution reproduces the exact result with a maximal 
deviation which is slightly above 1% as /3 2 — > 8ir. Entering the attractive region below the free 
fermion point this contribution falls down quite rapidly as the particles Aj become heavier than 
the lowest breathers (see (|4.3|l ). The inclusion of the first few breather states, however, gives a 
quite accurate result also in the repulsive region (Fig. 7). 

The same qualitative pattern can be observed in Figs. 8, 9 and 10 showing results for the 

rati ° S x approx 

= (6.9) 
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Figure 8: The ratio for the operators fij coming from Eqs. (|fi.4|) and (jdlj) . The squares 

indicate the contribution of the states AjAj only, the diamonds the inclusion of the states B2 
and B\B\, the stars the inclusion of the state -B4. 




2 4 6 8 Jl 



Figure 9: As in Fig. 8 for the operator V* . 
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Figure 10: As in Fig. 8 for the operator £ . 



between the approximations for the scaling dimensions obtained inserting a truncated spectral 
series in (|d4j) and the exact values. Quantitatively, the absence of the ultraviolet suppressing 
factor in (|6.4j) as compared to 1)6.3(1 obviously leads to poorer results for a given level of trun- 
cation. Of course, the accuracy can be improved by including more states in the spectral sum. 
Here we simply observe that the increasing deviation from the exact value for larger values of /3 2 
in the repulsive region is due to the increasing ultraviolet singularities of the exact correlators 
entering the sum rules, what makes increasingly important the contribution of the short dis- 
tances to the integrals. The result for = re is plotted up to the free fermion point because the 
corresponding integral in (J6.4|) diverges for j3 2 > 4ir (see |28j ) . The integral entering the com- 
putation of r-p* becomes divergent at f3 2 = 8tt, what helps understanding why the two-particle 
approximation is particularly poor as this point is approached. 



7 Universal ratios 

The ability to compute the correlation functions allows the evaluation of the canonical thermo- 
dynamic observables. We will consider the magnetisations 

M a = | fa) | (7.1) 
M P = \(V)\ (7.2) 
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the specific heat 
the susceptibilities 



C = / d 2 x (S (x)£(0)) c 



Xa= d 2 x (aj(x)aj(0)} ( 



Xl2 



XV 



d 2 x{a 1 (x)a 2 (0)) l 
d 2 x{V(x)V(0)) c 



the second moment correlation lengths 



-2nd 



1 



4Xa 
(J- 



d 2 x \x\ 2 (<jj(x)<jj(0)) c J 



1/2 



d 2 x\x\ 2 {V(x)V{0)) ( 



1/2 



and the exponential correlation lengths defined as 



kl 



lim (£r,-(a;)(T 3 -(0)} £ 



lim (V(x)V(0)) c ~ e 



-|*I/€p 



(7.3) 

(7.4) 
(7.5) 
(7.6) 

(7.7) 
(7.8) 

(7.9) 
(7.10) 



These quantities can be evaluated on the two sides of the critical line with continously varying 
exponents. In the disordered phase II we decompose the above correlators onto the form factors 
of section 5. Duality is exploited to obtain the results for the ferromagnetic phase I: the operators 
are replaced by their duals and the correlators are again decomposed over the form factors of 
section 5; the observables above do not depend on which of the four ferromagnetic ground states 
is selected by spontaneous symmetry breaking. 

For a given observable we denote by J- + its limit toward a given point on the critical line 
along a path in phase II, and by T~ the limit toward the same point along the dual path in phase 
I. Dimensionless numbers independent on metric factors can be obtained suitably combining the 
limits toward the same fixed point of different observables. These numbers are universal and 
characterise the scaling region around the critical line with continously varying exponents. 

Some of these universal combinations can be determined exactly. It follows from the spectral 
decomposition (|dl[) that the exponential correlation length is simply the inverse mass of the 
lightest asymptotic state coupling to the operator <£. In the disordered phase the lightest state 
is Aj for cTj-; for V we have instead A\A 2 in the repulsive region j3 2 > 4ir (i.e. £ > n) and B\ 
in the attractive region. In the ferromagnetic phase, both aj and V (or equivalently fij and V* 
in the disordered phase) couple to AjAj for (5 2 > 8n/3 (i.e. £ > 7r/2) and to B2 below this 
threshold. The following universal ratios then follow from (|4.3JI 



So- 




(7.11) 
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S<7 



1, f > 7T 



1/sinf, f<£<vr (7.12) 
2 cos | , £ < f 

2, f > 7T 



i ) am 3. 



(7.13) 

2 sin | , £ < 7r . 

The fact that the energy operator £ is odd under duality and the specific heat is bilinear in 
£ implies 

C + /C~ = 1 . (7.14) 

The combination 

C + (C) 2 = -(l-a)(2-a)^, (7.15) 

where a = 1 — -| is the specific heat critical exponent and (0) is given in 1)6.8(1 . is also exact. 

The universal ratios involving the susceptibilities and second moment correlation lengths 
cannot be computed exactly. In Table 4 we list the results provided by the form factor approach 
including in the spectral series all the one and two-particle states (two-particle approximation). 
The quantities R a and R-p are defined as 

An important indication about the size of the error involved in the two-particle approximation 
comes from the comparison with the exact results \A0\ avalaible for the point f3 2 = 47r, where 
the system reduces to two decoupled Ising models. We stress that, although the theory at this 
point is free, the opearators (Xj and V belong to the non-trivial sector and have non-zero form 
factors on an arbitrary number of particles. Hence, the results obtained for their correlators are 
representative of what happens at generic values of (5. Table 5 shows that the error of the two- 
particle approximation is extremely small (less than 0.1%) for the ratios involving er, only, while 
it grows to order 1% for the ratios involving V . This fact has a very clear origin. The spectral 
representation we use for the correlators is a large distance expansion and when we truncate it 
to obtain approximated results we make an error on the 'short' distances. Hence, the error on 
the integrals over all distance scales grows with the strength of the ultraviolet singularities of 
the correlators, namely with the scaling dimensions of the operators. The scaling dimension of 
V is twice that of Oj at f3 2 = 47r and this explains the two different error scales. 

On these grounds we expect that the error on the ratios involving only aj will stay quite 
small along the whole critical line, as a consequence of the fact that the scaling dimension of 
this operator does not depend on (3. Concerning the ratios involving V, the value X-p = -j^ 
suggests that the error will be of the same size of that of the a-ratios at 1 = 2ir and then will 
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0.9997 


3.002 


33.26 


0.4356 


1.931 


0.5438 


1.548 


5.175 


0.5724 


3.3 


0.9998 


3.046 


34.38 


0.3769 


1.948 


0.5224 


1.513 


4.847 


0.5323 


3.4 


0.9999 


3.081 


35.33 


0.3184 


1.963 


0.5026 


1.479 


4.538 


0.4950 


3.5 


0.9999 


3.110 


36.10 


0.2608 


1.975 


0.4840 


1.445 


4.247 


0.4602 


3.6 


1.000 


3.131 


36.71 


0.2045 


1.984 


0.4668 


1.412 


3.972 


0.4278 


3.7 


1.000 


3.146 


37.16 


0.1500 


1.991 


0.4506 


1.380 


3.711 


0.3976 


3.8 


1.000 


3.155 


37.47 


0.09763 


1.996 


0.4355 


1.349 


3.464 


0.3694 


3.9 


1.000 


3.161 


37.64 


0.04756 


1.999 


0.4214 


1.319 


3.228 


0.3430 


4.0 


1 


3.162 


37.70 





2 


0.4082 


1.291 


3 


0.3183 


4.1 


1.000 


3.161 


37.65 


0.04494 


1.999 


0.3959 


1.264 


2.782 


0.2952 


4.2 


1.000 


3.157 


37.50 


0.08721 


1.996 


0.3843 


1.238 


2.579 


0.2736 


4.3 


1.000 


3.152 


37.28 


0.1268 


1.992 


0.3735 


1.213 


2.389 


0.2533 


4.4 


1.000 


3.145 


36.98 


0.1637 


1.986 


0.3634 


1.189 


2.211 


0.2343 


4.5 


1.000 


3.138 


36.63 


0.1981 


1.979 


0.3539 


1.165 


2.044 


0.2165 


4.6 


1.000 


3.130 


36.23 


0.2299 


1.971 


0.3450 


1.142 


1.887 


0.1998 


4.7 


1.000 


3.122 


35.78 


0.2594 


1.961 


0.3366 


1.120 


1.740 


0.1841 


4.8 


1.000 


3.114 


35.30 


0.2866 


1.950 


0.3288 


1.099 


1.603 


0.1694 


4.9 


1.000 


3.105 


34.80 


0.3117 


1.937 


0.3215 


1.078 


1.473 


0.1556 


5.0 


1.000 


3.097 


34.28 


0.3348 


1.924 


0.3147 


1.059 


1.352 


0.1427 


5.1 


1.000 


3.090 


33.74 


0.3560 


1.909 


0.3082 


1.040 


1.238 


0.1306 


5.2 


1.000 


3.082 


33.18 


0.3754 


1.893 


0.3022 


1.022 


1.132 


0.1193 


5.3 


1.000 


3.075 


32.62 


0.3932 


1.877 


0.2966 


1.005 


1.032 


0.1087 


5.4 


1.000 


3.069 


32.05 


0.4095 


1.859 


0.2914 


0.9894 


0.9383 


0.09876 


5.5 


1.000 


3.063 


31.48 


0.4244 


1.840 


0.2865 


0.9742 


0.8509 


0.08950 


5.6 


1.000 


3.057 


30.91 


0.4380 


1.821 


0.2819 


0.9599 


0.7691 


0.08086 


5.7 


1.000 


3.052 


30.33 


0.4503 


1.800 


0.2776 


0.9463 


0.6929 


0.07280 


5.8 


1.000 


3.047 


29.76 


0.4616 


1.778 


0.2736 


0.9336 


0.6218 


0.06530 


5.9 


1.000 


3.043 


29.18 


0.4719 


1.756 


0.2699 


0.9216 


0.5557 


0.05834 


6.0 


1.000 


3.039 


28.60 


0.4812 


1.732 


0.2664 


0.9104 


0.4944 


0.05188 



Table 4: Two-particle approximation for the universal ratios along the Ashkin- Teller critical line. 
The relation between (3 and the values of the lattice couplings J and J\ at the corresponding 
critical point is provided by Eqs. (|3.15j) and Q2.3J) . The values /3 2 /vr = 2,8/3,4 correspond to 
the 4-state Potts, decoupled Ising and Fateev-Zamolodchikov models, respectively. 
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Ratio 


Two-particle 


Exact 




approximation 




e nd+ /e 


1 


0.999598087.. 




3.1623 


3.16249504.. 


Xcr / Xcr 


37.699 


37.6936520.. 




2 


2.00163051.. 


& nd+ /C 


0.408 


0.40656.. 


/s-p 


1.291 


1.3088.. 


Xp/Xp 


3 


2.9108.. 


#P 


0.318 


0.32104.. 



Table 5: Universal ratios at the Ising decoupling point (3 2 = 4n. 



increase with (3 to values that (from the results of the previous section) should not exceed 10% 
at P 2 = 6tt. 

For (3 2 = 2n the results we obtained should reproduce those for the 4-state Potts model. In 
order to check this point we label a = (01,02)) with ai j2 = ±1, the four states in which each 
site of the lattice can be. Then we build out of the Ashkin-Teller spin variables o\ and <r 2 the 
site variable a = (a±, cr 2 ) and introduce the traditional Potts spin variables 

= ^[aiCTi(x) + q 2 ct 2 (x) + aia 2 <7icr 2 (x)] (7-17) 

satisfying J2a s a = 0. If we denote by |0 7 ) the ground state that spontaneous symmetry breaking 
has selected in the ferromagnetic phase I, we will have 

(0 7 |s Q (x)s a (0)|0 7 ) = -^(0 7 |<ti(x)cti(0) + 0- 2 (x)0- 2 (O) + <Ticr 2 (ir)<7io- 2 (0) + 2aia 2 ai(x)a 2 (0) 
+ 2a 2 cri(x)(JiCT 2 (0) + 2aiCT 2 (x)(7icj 2 (0)|0 7 ) . (7.18) 

Along the Potts trajectory the internal symmetry is enhanced to invariance under global per- 
mutations of the four colours and one finds the expression 

(0(i,i) \s a (x)s a (0)\0( 1A )) p u s = 

■^<°(i,i) |3<ri (rc)cr! (0) + 2(a x a 2 + a 2 + ai)cri(x)<j 2 (0)|0( 1)1 ))p otts = 

1(0(1,1)13(71 (z)<7i(0) + 2(4«5 a;(1;1) - l)<7i(aO<7 2 (0)|0 (lil) > Potta (7.19) 

explicitely showing that only the cases = 7 and ct / 7 are distinguished in the correlator 
(0 7 |sq,(x)sq,(0)|0 7 ). In the disordered phase mixed correlators vanish by symmetry and we have 

3 

(s a (x)s a (0)) Potts = — {Gl{x)(Jl{ti)) Potts ■ (7.20) 
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In the Potts model we define 4 the longitudinal spontaneous magnetisation 



M = (0 a \s a \0 c 



(7.21) 



the high-temperature susceptibility 

X = J d 2 x(s a (x)s a (0)) c , 
the low-temperature longitudinal and transverse susceptibilities 

Xl = J d 2 x (0 a \s a (x)s a (0)\0 a ) c 
Xt = d 2 x (0 Q |s 7 (x)s 7 (0)|O a ) c 



a 7^7 



(7.22) 

(7.23) 
(7.24) 



and the second moment and exponential correlation lengths £ 2nd and £ computed from the 
correlator (s a (x)s a (0)) at high temperature and (0 a \s a (x)s a (0)\0 a ) at low temperature. 
The relations between these quantities and the Ashkin- Teller observables follow from Eqs. 1)7. 17|) . 
()7,19|) and Q7.2UJI . In particular one obtains 



X + /Xi 



Xt/x'i 



1 + 2xi 2 /x- 
— \x\2IXa 



~ 4.02 



/3 2 =2n 



~ 0.129 



(7.25) 



(7.26) 



J /3 2 =2tt 



l + 2 Xl2 /Xa 

These results agree 5 with those of Refs. E2] where the amplitude ratios for the g-state Potts 
model were computed. This non-trivial check eliminates the doubt raised in [32] about the 
result of |31j for the ratio x + Ix~l i n the (/-state Potts model. For q = 3 full agreement with the 
theoretical prediction was found in the lattice studies of Refs. |33) 134] . 



8 Conclusion 

We have computed the universal ratios along the Ashkin- Teller critical line with continously 
varying exponents using the field theoretical description of the scaling limit provided by the 
sine-Gordon model. As discussed in the Introduction, these results can be tested through nu- 
merical simulation or series expansions on the lattice model. Up to now, lattice results for the 
universal ratios along the Ashkin- Teller critical line exist only for the particular cases of the 
Ising decoupling point (where the exact values are known, see Table 5) and of the 4-state Potts 
model |351 134] . In the latter case the lattice results are compatible with the field theoretical 
predictions but are affected by large uncertainties originated by logarthmic corrections to scaling 
|361 137] coming from the marginal operator that is responsible for the end of critical line at the 
4 We drop the subscript Potts on the correlators below. 

5 There is a slight deviation from the values quoted in Refs. |31l I32| due to the fact that the contributions of 
the states AjB-z and B2B2 had been neglected in those works. 
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Potts point. Reducing the error bars at this point as well as obtaining estimates at other points 
along the Ashkin- Teller critical line appear as challanging tasks for future lattice studies. 
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